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We determine the mass of the charm quark (m c ) from lattice QCD with two flavors of dynamical quarks with a 
mass around the strange quark. We compare this to a determination in quenched QCD which has the same lattice 
spacing (0.1 fm). We investigate different formulations of the quark mass, based on the Vector Ward Identity, 
PCAC relation and the FNAL heavy quark formalism. Based on these preliminary results we find no effects due 
to sea quarks with a mass around strange. 



1. INTRODUCTION 

Quark masses are fundamental parameters of 
the Standard Model but due to confinement they 
cannot be measured directly by experiment and 
therefore any mass quoted always depends on the 
scheme (and scale) in which it was computed. 
The present world average value for m^ s (m c ) has 
an error of approximately 17% pQ. This should 
be compared to the value quoted for the experi- 
mental mass of the D meson which is accurate to 
within less than 0.1%. 

The mass of the charm quark is a difficult quan- 
tity to compute since it is light enough to be chal- 
lenging for heavy quark methods and yet heavy 
enough such that the lattice spacing may be too 
coarse for the Compton wavelength. There has 
recently been several calculations of this quantity 
from the lattice |2l.'il H5j . In particular, a detailed 
study by Rolf and Sint [3] has been presented in 
the continuum limit, thus eliminating the error 
arising from lattice artifacts. However, the ef- 
fects of including quark loops within these calcu- 
lations has so far been omitted. This work is an 
attempt to estimate these effects and obtain a re- 
sult in full QCD (albeit at fixed lattice spacing). 
In the interest of reducing errors resulting from 



a finite lattice spacing, we investigate several dif- 
ferent formulations of m c on the lattice. 

We have previously published results from this 
data set on the hyperfine and P-wave mass split- 
tings H . The mass splittings between the P-wave 
and S-wave mesons have been compared to the 
mass of the D*j(2317) + meson recently discov- 
ered by BaBar. 



2. DEFINING THE QUARK MASS ON 
THE LATTICE 

In this section we outline the definitions of bare 
quark mass used in this study and discuss some 
aspects of the renormalisation and improvement 
procedures. We use three different definitions, the 
first of which is the vector Ward Identity ( VWI) , 
given by 
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where k is the hopping parameter, the subscript 
h denotes the heavy quark and the subscript crit 
denotes the value of k which corresponds to zero 
quark mass. 

The second definition arises from the PCAC 
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relation (AVWI) and is given by 



(2) 



where Q and g label the heavy and light quark 
respectively, A is the local axial current and P is 
the pseudoscalar density. The axial current has 
been 0(a) improved. 

The scaling studies of Rolf and Sint |3] 
show that 0(a 2 ) lattice artifacts of the non- 
perturbative renormalisation scheme can be large. 
As some of the required coefficients have not 
been computed non-perturbatively at the pa- 
rameter values used in this study, we use the 
(boosted) perturbative values, thus the lattice ar- 
tifacts are 0(a s a). As the lattice spacing is rel- 
atively coarse, we utilise an alternative definition 
of m c . This is obtained by using the bare value 
obtained from the VWI, which is then identified 
as the "rest mass" in the FNAL formalism [7] . At 
tree level, this is defined as 



ami = log(l + amy)- 



(3) 



We use the one loop expression to connect the 
vector mass to the quark mass from Ref. In 
this work we use the rest mass of the hadron to 
compare with experiment, in the FNAL formal- 
ism the kinetic mass of the hadron is used. This 
is planned in future work. 

Matching onto the MS scheme is performed at 
one loop in perturbation theory for all definitions 
of the quark mass at the scale q = 1/a and the 
associated systematic uncertainty is estimated by 
matching at q = it /a. 

3. DETAILS OF THE COMPUTATION 

3.1. Lattice specifications 

There are four ensembles of gauge configura- 
tions available for computing the meson corre- 
lation functions but here we present details and 
results for the dynamical and coarsest quenched 
set only. The dynamical ensemble was gener- 
ated with two degenerate flavors of sea quarks 
with a mass around the strange quark. Further- 
more, the lattice spacing for the dynamical set 
was matched to that of the coarsest quenched en- 
semble (presented here) to facilitate the study of 



effects due to sea quarks without ambiguities aris- 
ing from different lattice spacing errors. All en- 
sembles were generated using a Wilson gauge ac- 
tion and the non-perturbatively O(a) improved 
fermion action was used to generate the quark 
propagators. These details are summarised in Ta- 
bleland further details of the procedures for gen- 
erating the dynamical ensemble and matching to 
the quenched data set can be found in Ref. UJ. 



Table 1 

Lattice parameters 



(/?; K sca) 


volume 


a- L (GcV) 


con figs 


(5.93,0) 
(5.2,0.1350) 


16 3 x 32 
16 3 x 32 


1.664$) 
1.716(44) 


347 
395 



The scale has been set using the "Method of 
Planes" technique which corresponds to setting 
ro = 0.55 fm as the central value. The asso- 
ciated systematic error is then estimated using 
ro = 0.5 fm. 

3.2. General approach 

Having computed the correlators, the meson 
mass is obtained as a function of quark mass as 
follows. We fit to a 2 by 2 matrix of correlators 
(elements of which are combinations of local and 
fuzzed operators) and the fit region is obtained by 
studying the effective mass plots and monitoring 
X 2 /dof. 

By extrapolating (interpolating) mii g ht to m u 
(m s ), the D and D s meson masses are obtained 
as functions of m-heavy We then plot the data 
using the following interpolation formula 



aM H = a + aiTOheavy 



(4) 



The experimental results for Md and Md s are 
then used to read off the bare value of m c . 

4. PRELIMINARY RESULTS 

Table[5]contains preliminary results for m c (m c ) 
in the MS scheme for the quenched and partially 
quenched sets with the same lattice spacing. The 
FNAL result presented here is obtained from the 
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rest mass definition and the errors quoted are sta- 
tistical and systematic respectively. These results 
have been converted into Renormalisation Group 
Invariant (RGI) masses and plotted against the 
lattice spacing in Fig. ^ with some previous re- 
sults shown for comparison. 

Table 2 

Preliminary results for m^ s (m c ). 



definition 


NF 


m c (m c ) (GeV) 


AVWI 


2 


1.705(5)(115) 


VWI 


2 


0.905(5)(130) 


FNAL 


2 


1.228(3)(120) 


AVWI 





1.606(4) (96) 


VWI 





0.864(5) (95) 


FNAL 





1.261(4) (95) 



AWI (NF=2) UKQCD 
VWI (NF=2) UKQCD 
FNAL (NF=2) UKQCD 
AWI (NF=0) UKQCD 
VWI (NF=0) UKQCD 
FNAL (NF=0) UKQCD 
VWI (Z A ) Rolf and Sint 
AWI Rolf and Sint 
VWI Rolf and Sint 
AWI Becirevic et al 
VWI Becirevic et al 
Juge 



(a/r ) 

Figure 1. Recent data for the charm mass in the 
RGI scheme as a function of the lattice spacing. 

5. CONCLUSIONS 

There is no clear evidence of the effects of 
the sea quarks on the mass of the charm quark. 
This is almost certainly due to the unphysically 
heavy sea quark masses. However, large system- 
atic uncertainties make it difficult to compare the 
quenched and dynamical results. One major hur- 
dle in reducing the systematic error is the size of 
the lattice spacing. The dynamical results have 
been obtained at a fixed lattice spacing which is 



coarse. In comparison with the quenched case, 
it is very expensive to reduce the spacing in un- 
quenched simulations and therefore we are unable 
to take the continuum limit which would enable 
us to control the systematic error. In response 
to this, we are investigating which formulation of 
m c on the lattice minimizes discretization errors. 

The mass independent renormalisation scheme 
suffers from large lattice artefacts of 0((am) 2 ) 
for quark masses at coarse lattice spacings [3]. 
Some of the coefficients required for this scheme 
have large O(a) ambiguities themselves. More- 
over, not all of these coefficients are known non- 
perturbatively for the dynamical case and so we 
have resorted to one-loop perturbation theory in 
this scheme. The only way to control errors fur- 
ther is to take the continuum limit which, as dis- 
cussed, is unfeasible with the current generation 
of computers. An alternative approach is that of 
heavy quark methods such as the FNAL formal- 
ism, which removes mass dependent lattice arte- 
facts of 0(a s (am) n ). With careful application of 
this formalism, we may be able to make a compar- 
ison between the quenched and dynamical case. 
This is the subject of future work, as well as tak- 
ing the continuum limit in the quenched case to 
control remaining lattice artefacts. 
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